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The BFKL fan diagram equation for the scattering on the nucleus is solved numer- 
ically with the eikonalized initial condition and for a realistic nuclear density. The 
gluon density has a soliton-like form in the log q — y space. Inclusive cross-sections 
for jet production in hA and AB collisions are calculated. 



1 Introduction 

In the colour dipole approach with N c — > oo 0! the interaction with a heavy 
nucleus is exactly described by a sum of fan. diagrams constructed of BFKL 
pomerons, each of them. splitting into two I3'D. The equation for this sum has 
been constructed in aaa. Let the forward scattering amplitude on the nucleus 
at fixed impact parameter b be 

A(y, b) = 2is J d 2 rp{r)<5>{y, b, r). (1) 

where p is the colour density of the projectile and $ represents the sum of 
all BFKL fan diagrams attached to it. Then function <f> — $/(27rr 2 ) in the 
momentum space satisfies a non-linear evolution equation in the rapidity 

(d/dy + H)cf ) (y, q ,b) = -q> 2 (y, q ,b), (2) 

where y = Y/a with a = a s N c /ir is a rescaled rapidity and H = In r 2 + In q 2 + 
const is the forward BFKL Hamiltonian. Eq. (2) should be solved with an 
initial condition which follows from the Glauber form of $ at zero rapidity: 

$(0, r, b) = 1 - exp ( - 8tt 2 AT(b) J d 2 r'G(0, r, r')p N (r'fj . (3) 

Here G is the forward BFKL Green function, T is the nucleus profile function 
and pn is the colour density of the nucleon from the target nucleus. Eq. (2) 
was studied perturbatively in u, by asymptotic estmates in □ and finally solved 
numerically in □ under the simplifying assumptions that the nucleus has a finite 
radius and that multiple Glauber interactions in (3) can be neglected, which 
is justified if A 1 / 3 « N c . 
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At small q the solution of (2) behaves like const — log q and consequently 
&{y, r, b) — > 1 at y — > oo and fixed r and b. This corresponds to the saturation 
at high rapidities of the cross-section for the scattering of a colour dipole on the 
nucleus to its unitary black disk limit. However this does not imply that the 
structure function of the nucleus saturates, since the number of colour xlipoles 
generated by the virtual photon is unlimited. In fact, as argued in vsu and 
confirmed by numerical calculations in it grows with y. The gluon density 
of the nucleus in the combined momentum-impact parameter space 

dxG(x,q,b)/d 2 qd 2 b = d 2 4>/{2ira{d\nq) 2 ) = h/(2%a), y = -\nx (4) 

was found to be a soliton wave in y— In q space moving towards higher rapidities 
with a constant velocity and preserving its nearly Gaussian shape. A fit to 
numerical data with £ = In q gives 

h(y > q ) b) = h e- a ^-^y^ 2 ) £o = -3.11 + (2/3) InB + A y, (5) 

where Ao = 2.3a, ho ~ 0.3 and a ~ 0.3 are universal. Dimensionless parameter 
B = na 2 AT{b)R 2 N where Rn is the nucleon radius. Actually the asymptotics 
of h at high |£| is not Gaussian but exponential: at £ — > oo h ~ exp(— 1.3£) 
and at £ — > — oo h ~ exp 2£. 

Here we first report on the improved numerical solution of Eq. (2) which 
takes into account multiple Glauber collisions in (3) and a realistic form of 
nucleus. Next we consider simplest of the production processes: the single 
inclusive jet production in hA and AB collisions and double inclusive jet pro- 
duction in hA collisions. 

2 Improved nucleus structure function and gluon densities 

Our first improvement consists in taking into account the full Glauber series 
in the initial value (3). Since at large b and consequently small T(b) the 
contribution of the multiple scattering is neglegible, we used the same simplified 
profile function corresponding to a finite nucleus as in It was noted in U 
that the solution of (2) very quickly forgets its initial form as y grows, so that 
one expects to notice the influence of the glauberization of the initial function 
only at the initial stage of the evolution. This is confirmed by numerical 
calculations. The gluon densities obtained from the full initial function (3) 
and from its single rescattering term only differ at small y and q, where the 
former tends to zero and the latter tends to a constant. This difference turns 
out to be completely washed out already at y = 1. For y > 1 the two gluon 
densities are identical for all values of q. So our first conclusion is that the 
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Glauber rescattering at small y has practically no influence on the behaviour 
of the interaction with the nucleus at high rapidities. 

Next we studied the contribution of very peripheric collisions for realistic 
nuclei with exponentially falling profile function. The interest in this effect is 
due to the fact that this contribution is not damped by the non-linear term 
and grows exponentially with y until (f> ~ 1. If T(b) ~ exp(— b/ Ra) then this 
happens at b ~ 60 = Ra^V where A = 4 In 2 is the BFKL intercept. The 
total cross-section for the scattering of a colour dipole on the nucleus at b > bo 
and y >> 1 then grows linearly with y. Corespondingly one expects that the 
structure function will grow with y faster than linearly. These expectations are 
also confirmed by our numerical calculations. At large values of 1/x — expY 
(~ 10 15 -j- I0 20 ) the found structure functions of Pb with a constant density 
inside a finite sphere (CD) and with the Woods-Saxon nuclear density (WS) 
can be fitted by 

Ff D ~ Q 2 (332y - 73 In Q 2 - 893), 

F^ s ~ Q 2 (12 Ay - 2.58 In Q 2 - 29.0) 2 , 

where y = aY and Q 2 is in (GeV/c) 2 . Thus, to a good approximation, the 
structure function of a realistic nucleus grows with 1/x as In 2 (1/x). 

Finally a few words about our definition of the gluoii density in relation to 
the saturation phenomenon as discussed by A.Mueller □. We define the gluon 
density as essentially the quantity to be integrated with the quark loop to ob- 
tain the virtual photon scattering cross-section and structure function in the 
low x kinematical region. This definition is appropriate for this kinematical 
region, without strong ordering in the transverse momentum, and cannot be 
directly applied to the large a; region where such ordering takes place. This 
circumstance has to be taken into account when comparing our gluon density 

fith differently defined ones. In particular the gluon density introduced in 
via the interaction of a "gluonic current" with a nucleus is a totally dif- 
ferent quantity, more appropriate for the large x region. One can show that 
A.Mueller's density can be expressed in terms of our density (4) and that its 
gturation properties obtained in □ follow from the found solitonic form of (4) 

3 Single jet inclusive production 

Single jet production inclusive cross-sections in hA collisions are obtained from 
the imaginary part of forward scattering amplitude A (Eq. (I)) by a substitu- 
tion of one of the BFKL Green functions in the fan diagrams 

G(Y) —> G(Y - y)Vk(r)G(y), Vu(r) — (4N c a s /k 2 ) A e lkr A (6) 
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(the arrows shows the direction of differentiation). Due to the AGK rules con- 
tribution of all such substitutions below the uppermost spliiting point cancel 
so that we are left with 

da/dyd 2 kd 2 b = 2(pG(Y - y)V k $(y)) = (16ir 2 a/k 2 )(p[G(Y - y) A]e lkr h(y)), 

(7) 

where p and $ are from Eq. (1), (...) means integrating over the gluon relative 
transverse coordinates or momenta and the relation between $ and <fr and the 
explicit form of Vk have been used to obtain the second equality. 

Using the form of h obtained numerically we find jet multiplities at large 
y and Y — y 

p A (Y, y, k) = (c/k 2 )A 2 / s e AY -*y/Vy(Y-y) (8) 

where A = 4 In 2, e= 0.39, c is a known numerical constant and all rapidities 
are scaled with a. The inclusive cross-section integrated over k evidently di- 
verges at low k 2 . Physically relevant results correspond to jets with not too 
small transverse momentum k > k m i n . 

The A dependence is different both from the eikonal and local pomeron fan 
diagrams predictions. The k dependence is 1/k 2 up to lnfc ~ Ay in contrast 
to the hh scattering where it is cut by the damping factor exp(— cm 2 k/y). In 
this way we observe the " Cronin effect" : the distribution for the nucleus target 
is flatter than for the hadronic one. 

These.jaesults are easily generalized to nucleus- nucleus scattering. As was 
shown in U, for them the single inclusive cross-section is given by the sum of 
fan diagrams connecting the produced jet with the both participant nuclei. So 
at fixed impact parameters of the nuclei b A and b B the inclusive cross-section 
is given by an obvious generalization of Eq. (7): 

da/dyd 2 kd 2 b A d 2 b B = 2(<f> A (Y-y)V k $(y)) = (32ir 3 a/k 2 )(h A (Y-y)e lkr h B (y)) 

(9) 

Using our results for the gluon densities, for lnfc << £o, y ~ Y/2 and 
identical nuclei, we obtain rough estimates ~ (l/k 2 )A 4 ^ 9 e AlY for the inclusive 
cross-section at a given k and ~ ^4 2 / 3 e t A i y for the cross-section integrated 
over fc > k m in with Ai = 2.4. Note that from the latter estimate it follows 
that the multiplicity density at a given y behaves as A 4 / 3 , that is, similarly to 
the eikonal result. 

These estimates are supported by our numerical calculations. For Pb-Pb 
collisons at center rapidity with a = 0.16 we obtained values for the multiplicity 
which are well described by a fit 

M (0) = 0.092 e °- 59y (l + 35. e -° 15F ). 
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Note that with a = 0.16 the BFKL pomeron intercept is A = 0.42, so that the 
multplicities grow faster than the BFKL pomeron. The A dependence of /i(0) 
at a given Y is almost perfectly represented by a power factor /x(0) ~ A@( Y \ 
the power 0(Y) slowly rising with Y from 1.25 at Y = 10 -f- 20 to 1.29 at 
Y = 40. So it seems that the difference in the A dependence with the eikonal 
prediction {(3 = 4/3) is gradually disappearing at high enough rapidities. 

The form of the distribution in y turns out to be practically independent 
of the atomic number. With the growth of energy the distribution becomes 
more strongly peaked at the center. 

4 Double inclusive cross-sections 

For the nucleus-nucleus collisions the double inclusive cross-sections cannot 
be expressed via sums of fan diagrams for the participant nuclei, but require 
summation of diagrams of a more complicated structure c3. For this rea- 
son we restrict ourselves to the case of hA scattering. The double inclusive 
cross-section is obtained from the imaginary part of the amplitude (1) by two 
substitutions (6) for the two produced jets with rapidities yi$ and transverse 
momenta fci^- The AGK rules leave two contributions corresponding to the 
two substitutions either in the uppermost pomeron or in the two pomerons im- 
mediately after the first splitting. These two contributions rise with the overall 
rapidity Y as e AY and e 2AY respectively, so that the second one is dominant 
at high Y. 

Study of this contribution shows that the multiplicity ratio 



R= n(yi,ki;y2,k 2 ) = c (Y - yi)(Y ~ y 2 ) 
fx(y!,ki)fj,(k2,y2) \j 2Y-y 1 -y 2 ' 

where c is a known constant proportional to a s , does not depend on fci^, grows 
with Y at center rapidity and stays constant in the diffractive regions. With 
a s = 0.2 at y = Y/2 we find R = 0.19y/Y so that correlations change from 
negative to positive at Y ~ 25. 

5 Conclusions 

Improved numerical calculations confirm the validity of approximations made 
in solving Eq. (2) in B. The nuclear structure function seems to grow with 
rapidity also for a realistic nucleus with a profile function exponentially falling 
at large impact parameters, the growth being faster than for a finite nucleus. 

The inclusive jet production cross-sections grow rapidly with energy as the 
pomeron itself or even faster (for AA collisions). Their A- dependence results 
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close to the eikonal predictions in evident contrast to the earlier ones made for 
fan diagrams with a local pomeron. 
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